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^ ■ Abstract 

■ Nematodynamics is the orientation dynamics of flowless liquid-crystals. We 
^ ■ show how Euler-Poincare reduction produces a unifying framework for various 

^ . theories, including Ericksen-Leslie, Luhiller-Rey, and Eringen's micropolar theory. 

" I In particular, we show that these theories are all compatible with each other and 

■ some of them allow for more general configurations involving a non vanishing 
disclination density. All results are also extended to flowing liquid crystals. 
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1 Introduction 



1.1 Ericksen-Leslie theory vs. Eringen's micropolar theory 

A flowless nematic liquid crystal possesses orientational dynamics that involves the di- 
rector variable n(x, t) along with its associated angular velocity. The director n is a unit 
vector, i.e., |n| = 1, that is unsigned, so that n is identified with — n (the head is iden- 
tified with the tail). This orientational dynamics is often called nematodynamics and 
its equations are subject of continuing research. Depending on whether inertial effects 
are considered or not, different theories are obtained. When these effects are neglected, 
dissipation plays a fundamental role and the problem is often approached from a kinetic 
theory point of view [3]. Although most of realistic situations involve negligible inertia, 
nematodynamics theories are best understood when inertial effects are taken into ac- 
count. For the purposes of this paper, we shall consider the ideal cases when dissipation 
is absent and inertial effects are not negligible. For example, in the ideal case of zero 
dissipation, the celebrated Ericksen-Leslie equation 

has been given a well defined Hamiltonian structure; see [1]. In the above equation, J 
is the microinertia constant and the molecular field 

5F d 5F 



5n dx^ 6{dx^n) 

is expressed in terms of the Frank F energy, which in turn is given by 

11 1 
F(n, Vn) := (n ■ curln)^ +-J'rii (divn)^^+-J'r22 (n ■ curln)^^ + -A'33 ||n x curlnf^ , 

chirality splay twist bend 



where each term possesses a precise physical meaning, as indicated above. 

The Ericksen-Leslie equation for the dynamics of nematic liquid crystals is a widely 
accepted model since it has been experimentally validated through various measurements 
[2, 1, 22]. However, orientational defects (disclinations) , make this model unreliable. For 
example, if defects are present (singularities in the orientational order), uniaxiality may 
be lost and the model becomes biaxial thereby eliminating the usefulness of the director 
field n: one needs some other order parameter. 

There are other approaches that take into account defect dynamics. The most gen- 
eral one was proposed by Eringen [7] in the context of microfiuid motion that includes 
liquid crystals. The Eringen model admits molecular shape changes by introducing a 
microinertia tensor j, whose dynamics is coupled to the wryness tensor 7, which is often 
expressed in terms of (Vn) x n when disclinations are absent [7]. In this paper, the 
notation Vn stands for the vector valued one-form that is given by the derivative of the 
director n, viewed as a map with values in 5*^ C M^. 

Nematic liquid crystals are well known to be a typical example of microfiuids. Several 
unsuccessful attempts have been made to show how the Ericksen-Leslie (EL) description 
arises from Eringen's micropolar theory. The most notable one is the relation 7 = 
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(Vn) X n proposed by Eringen [7, formula (11.2)]; it does not yield the EL equations 
[18] as shown in [9, Theorem 8.11] by two different methods (symmetry considerations 
and a direct computation). No good formula exists that expresses 7 in terms of the 
director n. This difficulty arises even in the simplest ideal case when dissipation is 
totally absent. 

1.2 The gauge-theory approach 

Reduction theory has recently contributed to the understanding of defect dynamics 
[13, 9]; it leads directly to the gauge-theory approach [1] and applies to very general 
systems since it incorporates defect dynamics in various models, such as frustrated spin 
glasses [11, 1], for example. In this gauge theoretical setting, one interprets the wryness 
tensor 7 as the magnetic potential of a Yang-Mills field (i.e., a connection one-form) 
taking values in the Lie algebra so (3) of antisymmetric 3x3 matrices (usually identified 
with vectors in M^) of the proper rotation group S0{3). The one- form 7 is also known 
as 'spatial rotational strain' [13]; it usually expresses the amount by which a specified 
director field rotates under an infinitesimal displacement. The gauge potential 7 can be 
written in terms of a chosen basis {e^ | a = 1, 2, 3} of ~ so(3) as 

7 = Jidx' = -fiBadx'. 

Then, its corresponding magnetic vector field (the curvature of 7) is given component- 
wise by 

^* = e^''(5.7fc + 7, X7fc); (1-1) 

summation over repeated indices is understood and two-forms on physical space 
have been identified with vector fields. The absence of disclinations in the gauge theory 
approach (see [1]) is characterized by the vanishing of the magnetic field B and not by 
the vanishing of the vector potential 7. Therefore, if 7 7^ must be compatible with EL 
dynamics, then one requires that B = 0. In [ ] it was shown that B = comes down to 
the homogeneous initial condition Bq = (that is, for example, 7o = 0). If -B 7^ 0, then 
the gauge theoretical model confers the director formulation the possibility of including 
non-trivial disclinations. 

Although Eringen's micropolar theory appears to possess a gauge-theoretical for- 
mulation due to the gauge invariant relation Vn = n x 7, Eringen's final choice of 
identifying the wryness tensor 7 with (Vn) x n leads to problems arising from the fact 
that the latter expression does not transform as a magnetic potential under gauge trans- 
formations (i.e., it is not gauge- invariant; see [9, Lemma 8.10]). Still, Eringen's theory 
exhibits many properties analogous to gauge-theoretical models. In addition, the simul- 
taneous presence of the microinertia and wryness tensors in Eringen's model accounts 
for the interaction of molecular shape with a non-vanishing disclination density B. 

1.3 The role of reduction by symmetry 

The comments above are the main motivation for the present work. We shall employ 
Euler-Poincare variational methods to produce a framework that incorporates defect 
dynamics in continuum systems with broken internal symmetry (e.g., liquid crystals) 
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and show that Eringen's micropolar theory contains Ericksen-Leshe dynamics. The 
main idea is that the gradient of the identity 

n(x,t) = 63, 

relating director dynamics to the dynamics of the rotation matrix x(x, t) G 5*0(3) in 
EL theory, is 

Vn= (Vx)e3 = (Vx)x~'n, 

where 63 := (0,0, 1), the unit vector on the positive vertical axis. Next, note that the 
new variable 

7 = -(Vx)x-' (1-2) 
is a so(3)-valued connection one-form [13, 9]. 

Remark 1.1 (The hat map) In this paper we shall often make use of the following 
'hat map' isomorphism of Lie algebras ^: — )■ so(3) (see e.g., [20]) 

- _ _ 1 - 

holding for any vector 1/ G M'^ and its associated antisymmetric matrix D G so (3), so 
that uu = u X u for any vector u G M"^. In the formulas above, summation of repeated 
indexes is used. In particular, applying the inverse of the hat map to the relation (1.2) 
yields Eringen's wryness tensor 7jdx* [7] (see discussions below). 

The two relations above yield hence Vn = n x 7 which, when used in the EL equations 
introduce 7 as a new dynamical variable. Note that the choice 7 = ( Vn) x n satisfies the 
relation Vn = n x 7; however this latter relation determines 7 only up to a component 
parallel to n. Although this was noticed by Eringen in [6], he did not observe that, as 
a consequence, 7 is not a function of n alone; in fact all three columns of the matrix 
x(x, t) are needed to specify 7. 

The second main observation is that a different symmetry reduction of the same 
material Lagrangian yields a new set of equations for nematodynamics. This is due to 
a symmetry that seems not to have been exploited before. We will show that these new 
reduced equations are completely equivalent to the Ericksen-Leslie equations. However, 
this new system permits the description of disclination density dynamics, something 
that the Ericksen-Leslie equations could not handle explicitly. 

All the above considerations hold independent of the background fluid motion (the 
macromotion); they affect only the (orientational) micromotion. Because of this, we will 
treat for a while only liquid crystals with zero background fluid motion since this allows 
us to concentrate on the main points of the relation between the Ericksen-Leslie and 
Eringen models. At the end of the paper we shall quickly present liquid crystals with 
macromotion. 

2 Euler-Poincare reduction for nematic systems 

Nematic systems are continuum media in which each particle of the system carries an 
orientation, that is, a (time-dependent) rotational matrix x(x, t) G 5*0(3) is attached 
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at each point x in physical space M^. More specifically, we shall consider a medium 
occupying a bounded domain V C M^, so that x G P and, at each time t, x is a 
map from V to 50(3); the set of all such maps will be denoted by F{V, 50(3)). When 
dissipation is neglected, a nematic system can always be described by the Euler-Lagrange 
equations 

d5C _5C _^ 

(the notation x stands for the partial time derivative dtx) arising from Hamilton's prin- 
ciple 

6 C{x,x)dt = 0, 
Jti 

for any variation Sx(t) G J^(V, S0{3)) satisfying Sx(ti) = Sxih) = 0, where the La- 
grangian functional C involves a Lagrangian density =Sf as follows: 

£(X,X) = / ^(x(x,t),x(x,t))d=^x. 
Jv 

As usual, we think of x as a J^{T>, 5*0 (3))- valued curve, where P C M'^ is the physical 
domain of the liquid crystal and J^(T>, S0{3)) denotes the group (relative to point- 
wise multiplication) of smooth S'0(3)-valued functions on V. Thus the Lagrangian C is 
defined on the state space TJ^{V, 50(3)) and we denote a tangent vector to J^(V, 50(3)) 
atxby (x,x) eT^J-(P,50(3)). 

2.1 Continuum systems with rotational symmetry 

If Hamilton's principle were rotational invariant under right multiplication, that is 
^(X,X) = ^XX~^,XX~^) = >C(I,xx"^) for any x e J'{V,S0{3)), then one could define 
the angular velocity matrix u{'x,t) := x(x, t)x~"'^(x, t) G so(3) and the reduced La- 
grangian := C(J,xx^^) ill order to obtain the reduced Euler-Poincare variational 
principle 

S [\{i>)d^^ = (2.1) 

Jti 

for all variations Su = dt {{5x)x~^) + [(^x)x~\'^]- 

The variable v G J-'(P,M'^) is of paramount importance, first, because it is the only 
variable appearing in the reduced Lagrangian I. Second, the physical importance of the 
vector V is that it measures the angular velocity — (l/2)ejjfcXiaXaA; in the lab [spatial] 
frame and thus it is a physical observable. 

2.2 Nematic systems and the director 

The hypothesis of rotational invariance leading to the Lagrangian liy) is too restrictive 
in the context of nematic systems. This is due to the existence of a preferred direction 
of each particle encoded by the director variable n(x, t) = x(x, t)no(x), where no(x) is 
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the initial director at t = 0. The corresponding Lagrangian x) is usually written in 
terms of a free energy F, such that for nematic systems we have 

^X,x) = i [ |x(x,t)no(x)|'d3x- / F(x(x,t)no(x),V(x(x,t)no(x)))d=^x, (2.2) 
Jv Jv 

where J is the microinertia constant. Notice that this C is not ^-'(P, S'0(3))-invariant, 
that is, £(X) x) 7^ 'C(I) XX~^)- However, if we consider uq as an extra variable, we obtain 
the following invariance property 

^no{X,X) = >Cxno(I,XX~^) • 

Thus, if we define n := xno and £(t/, n) := £^no(I) the reduced variational 

principle reads 

6 i(iy,n)dt= { — ■5u + — -Sn]dt = 0, (2.3) 
Jh Jh V^*^ 6n J 

subject to the constrained variations that are obtained by a direct computation from 
the definition of the variables n and u, namely, 6n = {6x)x~^''^ ^ind 6u = dt {{5x)x~^) + 
[{Sx)x~^i I/]. This yields the reduced Euler-Poincare equations 

d 6i 6i Si 

at ov ou on (o a') 

dn ^ 
—— + n X 1/ = 0. 

dt 

Notice that the system is fully S'O (3)-invariant if and only if n = 0; otherwise one says 
the S'O (3) symmetry is broken by the presence of the special direction hq. 

The above equations provide the convenient setting for the derivation of the Ericksen- 
Leslie equation by taking the Lagrangian 

£(i/,n) = - I |i/xn|^d^x- I F(n,Vn)d3x, (2.5) 
2 Jv Jv 

computing its variational derivatives 

51 Si 

— = —Jn X (nx u) = Jnx dtn , — = —Jv x (i/ x n) — h = —Jv x dtu — h, 
ou on 

and applying nx to the z^-equation after using the Jacobi identity. 



2.3 Eringen's wryness tensor 

As we have seen, the free energy contains the derivatives of n. This is not a problem 
because Vn depends on the director n, so that the explicit dependence of F on Vn is 
irrelevant. For example, the molecular field is expressed as 

^_SF d 6F 
6n dx^ 6{d.j.in) ' 
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However, we want to emphasize that the spatial derivatives of n can be used to 
introduce an extra variable, as follows. Since 

Vn = V(xno) = (Vx)x"^n + xVno , 

it follows that if Uq is spatially constant, so Vng = 0, one can use the so (3)- valued one 
form 7 := — (Vx)x~^ on "D C to write diii = n x 7^, where the components of 7 
are 7^ = — (l/2)e"^'^7j'"^. If Vno 7^ 0, then there exist three non vanishing vectors 7oj 
such that diiiQ = no x 7Qj. This can be seen by recalling that the condition |no(x)p = 1 
yields diiiQ ■ no = 0, so that the three vectors 

7oi = diUo X no + A^no 

are defined up to an additive vector parallel to no- Then, from the previous equation, 
we obtain 

Vn = V(xno) = ((Vx)x"^ - X7oX~^) n =: n x 7 
where we have now defined 

7 := -(Vx)x"' + X7oX~' = x(Vx"') + X%X^'- (2-6) 

Following Eringen's work, we call the tensor 7 with components 7° the wryness tensor. 
It is important to notice that we sharply diverge here with Eringen's considerations, 
who defined 7 = (Vn) x n, a quantity that is not gauge invariant and would require 
7j ■ n = at all times. As we shall see below, it is precisely our definition of 7 that will 
lead to the main result: Eringen's micropolar theory contains EL nematodynamics as a 
special case. 

Returning to our definition (2.6) of 7 in terms of x, we notice that the initial condition 
7o appears precisely as a gauge transformation by x applied to the (Yang-Mills) magnetic 
potential 7o. This fact leads us to consider a nematic Lagrangian C{x, x) of the type 

^(x,x) = ^^lx(x,t)no(x)|^ d^x 

- / F(x(x, t)no(x), (Vx(x, t) - x(x, t)7o(x))no(x)) d^x . (2.7) 
Jv 

Again, the 5'0(3)-invariance is lost by the above total Lagrangian, that is £(x, x) 7^ 
C(I, XX~^)- However, if we consider no and 7o as extra variables, we obtain the following 
invariance property 

■^00,70 (X; '^^ -^11,7(17 XX ) 5 

where 

(n, 7) = (xno, x(Vx"') + X%X'') ■ (2.8) 

These definitions define an action of the gauge group J^{T>, S0{3)) on the Cartesian 
product T{V,S^) X n^{V,so{3)), where T{V,S^) is the space of director fields while 
Q^{V,so{3)) denotes the space Q^{V) ® 5o{3) (containing 7) of so(3)-valued one forms 
on V. See e.g., [20] for basic information on group actions. In this particular case, the 
X-action xno on the director no is linear, i.e., it is a representation. On the other hand. 
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the x-action xC^X^^) +X7oX~^ on the magnetic potential 7o is an ajfine action; see [9]. 
Notice that these group actions emerged naturally from the problem and the concept of 
a group action was not required at the beginning of our discussion, which started with 
a Lagrangian depending on the parameters no and 7o- 

Proceeding as in the previous case, we introduce the reduced Lagrangian £(y, n, 7) := 
Cn,^(J,XX~^) arising from the primitive Lagrangian £no,7o(X)X) i^i (2-7), and get the 
reduced variational principle 



6 i{iy,n,j)d^x = I — ■ 6u + — ■ 5n + — ■ SjA = 



(2.9) 



subject to the constrained variations 



{6n, 6j) = {{6x)x~'n, -V - [7, mx''] ) , 

6u = dt {{6x)x~') + [mx-\^] , 

obtained by a direct computation from the definition of the variables u = XX~^ 1 ^1 
and 7 given in (2.8). The variational principle (2.9) yields the reduced Euler-Poincare 
equations 

d 51 51 d 5i2 51 51 

— — = X h 7— h 7i X h n X — 

dt5v 5v ox' 57j 57 j 5w 

dtn + n X 1/ = 
dtli + 7i X 1^ + diU = 0, 

which are equivalent to the Euler-Lagrange equations for the primitive Lagrangian 
^nonoix^x)- These equations provide an alternative description of nematic systems. 
When 7q = and Vng = 0, the above system is fully equivalent to that obtained in 
the previous section. This follows from the fact that the Lagrangian £{x, x) material 
representation has never been modified and the only change was introduced through the 
initial condition Vng = no x 7q, which led to the introduction of 7. For example, upon 
replacing Vn = n x 7 (with 70 = 0) in the Frank energy we obtain a free energy $(n, 7) 
and the variational derivatives 

5$ 5F (5$ 5F 5F 



to produce 



5ji 5{din) ' (5n (5n * 5{din) 



^ 5$ 5$ (5$ 

n X h = di- h 7,' X t h n x — 

57, 57, 5n 



which leads to the compatibility between the two approaches. 

Remark 2.1 (Eringen's formulation of the wryness tensor) The explicit form of 
the wryness tensor lies at the origin of many problems that were encountered to show 
compatibility of Ericksen-Leslie theory with Eringen's micropolar theory. In particular, 
equation (2.6) was well known to Eringen (at least in the case when 7q = 0), who first 
defined the wryness tensor as 7 = —{'Vx)x~^- [6]; he noticed that 7 cannot be solved 
uniquely from the equation Vn = n x 7 and the identification 7 = Vn x n requires 
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7 • n = 0. However, in the same paper, it is claimed that the non- uniqueness of 7 
does not affect the expression of the free energy $, so that one can still use the relation 
7 = Vn X n. This statement is contradicted by the fact that the condition 7 ■ n = 
is not preserved in time and thus the quantity Vn x n cannot be identified with the 
dynamical variable 7 at all times. This point has been at the origin of all problems 
concerning Eringen's micropolar theory; problems that are now solved in the present 
paper, which shows how the micropolar theory extends Ericksen-Leslie nematodynamics 
by allowing for an arbitrary initial condition 7q. 

However, when 79 7^ 0, the relation between Ericksen-Leslie dynamics and Eringen's 
micropolar model is more subtle. Notice that the micropolar dynamics yields 

(^■^ - J^x^ (Vn - n X 7) = 

where n = xno and 7 = — (Vx)x^^ + XJoX"^ ^^^^ independent dynamical variables. In 
turn, this implies 

Vn - n X 7 = X (Vno - no x 79) . 

Then, if the initial conditions satisfy Vng = no x 7q, this relation is preserved by the 
dynamics and Eringen's micropolar theory is still equivalent to Ericksen-Leslie nemato- 
dynamics. On the other hand, when the initial conditions of micropolar dynamics are 
such that Vno 7^ no x 7q, then the two theories are not compatible. Notice that, when 
7o 7^ 0, a non vanishing magnetic Yang-Mills field = e*-''^(9j7ofc + 7oj x lok) defines 
the (initial) disclination density in the gauge theory of disclinations [4]. If the initial 
conditions are such that Bq = (e.g., 7oj = 0), then the disclination density vanishes 
at all times [9]. 



3 Broken symmetry and Noether's theorem 

3.1 Residual symmetries and isotropy subgroups 

In this section, we explain how the various ways in which symmetry can be used lead to 
different conserved quantities arising from Noether's theorem. These conserved quan- 
tities arise from the fact that, although the full 5*0(3) symmetry is broken in nematic 
systems (i.e., £(x, x) ^i^^XX"^)): residual symmetries are still possible. This means 
that there exists a Lie subgroup Go of J^{T>, S0{3) such that £(x, x) = '^{xXo^ ^ XXo^) 
for any xo ^ Go- This Lie subgroup leaving C invariant is known as isotropy subgroup 
and it is associated to the variable breaking the symmetry, e.g., the director no. In 
particular, the isotropy subgroup of no is defined as 

Gno = {Xo e J^{V, S0{3)) I xono = no} . 

One verifies immediately that the Lagrangian (2.2) is invariant under Gno- ^or simplicity, 
here we consider the case in which no = 63 so that Vno = 0. Then, 

G,, = J^{V,S0{2)), 
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the gauge group of planar rotations, Notice that, strictly speaking, one should make 
use of the (dihedral) group J^{'D, 5*0(2)) x Z2 to account for the equivalence no ~ —no; 
however, this reflection symmetry is left as an a posteriori verification in this treatment. 
This is consistent with most of liquid crystal theories. In principle, the discrete Z2- 
symmetry could also be taken into account by using sophisticated methods in reduction 
theory [21]. However, the use of these methods would require a much more involved 
treatment, which is left for future work. 

By proceeding analogously, one can also define the isotropy subgroup of 7o, upon 
using the affine x-action introduced previously, to get 



that is, the Lie group of spatially constant rotation matrices Xo ^ 5*0 (3). In more 
generality, the isotropy subgroup of the couple (63,70 — 0) given by Ge^ H G^y^^o = 
J^(T>, SO{2))nSO{3) = 5*0(2), i.e., the Lie group of spatially constant planar rotations. 

3.2 Geometric nature of conserved quantities 

Noether's theorem ensures that there are conserved quantities corresponding to each of 
the previous symmetry subgroups. The reduced equations provide a convenient setting 
for expressing these conserved quantities. In the simplest case of a fully J^iT), 5*0(3))- 
invariant system, one has no = (i.e., n = xno = 0) and the associated conserved 
quantity is x~^((5£/(5t'), that is. 



where i{u) is obtained from the Lagrangian appearing in (2.1) by using the hat 
map. The above relation is the specialization of a well known result in the geometric 
mechanics of invariant systems [20]. Notice that the momentum variable 6i/6u belongs 
to so(3)*, the dual space of the Lie algebra so(3) ^ M.^, and so does x'^^^/^^ (since 
the latter identifies the group coadjoint action Ad^-i di/dv). Thus one again identifies 
5o(3)* ~ 



As a general principle in geometric mechanics, conserved quantities belong to the 
dual space of the Lie algebra associated to the symmetry group leaving the system 
invariant. 



This means that if Hamilton's principle is invariant under the action of a Lie group G, 
then Noether's conserved quantity belongs to the dual space g* of the Lie algebra g of 



= {xo e HT^, 50(3)) I Xo7oXo ' + Xo(Vxo~') = 7o} ■ 



In the spatial case when 7o = (i.e., Vno = 0), its isotropy subgroup is 



G^„=o = S0{?>) 




G. 
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At this point we wonder what happens when the rotational symmetry is broken (i.e., 
n = 7^ 0) and the above conservation transforms into the relation 



d / 1 5£ \ „i / d 5^ 6i 



-1 



which arises from the first equation in (6.2). 

Following the general principle just stated, we seek a conserved quantity belong- 
ing to the dual Lie algebra of the isotropy group Gno, which leaves (2.2) invariant. 
Upon specializing to the case no = 63, we seek a quantity in the dual space of Qea = 
T{T>,so{2)) = C°°{V). A geometrically natural method to construct such quantity is to 
use the dual of the Lie algebra inclusion i{r) = (0,0, r) G J^(T>,M.^) ^ J-'(D,so(3)) for 
an arbitrary scalar function r G C°°(T>) ~ J^(T>,5o{2)). The dual operator i* is defined 
as f^fJ: ■ i{r) d'^x =: f^ri*{iJ,) d^x for an arbitrary vector function /x G J^(V,5o{3))* ~ 
^-"(1), M'^). Thus, ■i*(/x) = yU3, the third component of /x, and we obtain 

X TT: = U 17.) =X3^T- = n- — , 



where we have used Einstein's summation convention and the fact that ri; = Xik^^ok = 
Xik^ks = Xi3- III conclusion, the reduced Hamilton's principle (6.1) is accompanied by 
the conservation law 

d ^ S£\ 

n-— =0, (3.1) 



dt V Siy^ 

as can be easily verified from the Euler-Poincare equations (6.2). 

When the wryness tensor is also considered as a dynamical variable, then one seeks 
a conserved quantity belonging to the dual Lie algebra of the isotropy Lie subgroup 
Ge3,fg=o = 5*0(2). Using the same method as above, one considers the Lie algebra 
inclusion i{r) = (0, 0, r) G J^iT), M^) ^ J^(T>, 5o{3)) for any real number r G so(2). Then, 
the corresponding dual operator is = /p/i3(x)d^x and the conserved quantity 

associated to the reduced Hamilton's principle (2.9) is n ■ S£/Svd^:s. (see [10]), i.e., 

A/n.|^d^x = 0, (3.2) 
dt ou 

as long as no = 63 and 7o = 0. 



3.3 Conserved quantities are momentum maps 

In geometric mechanics, Noether's conservation laws are reformulated in terms of mo- 
mentum maps. These maps take values in the dual of a Lie algebra and always identify 
the conserved quantities of symmetric systems, i.e., systems whose Hamilton's principle 
is invariant under the action of a Lie group Gq. The various (and important) properties 
of momentum maps assume some basic knowledge of Lie group theory, which can be 
found in geometric mechanics books (e.g., [20]). One of the most important results in 
mechanics linking symmetry to conservation laws states that 

momentum maps always identify Noether's conserved quantities when they arise from 
the same symmetry group that leaves Hamilton's principle invariant. 
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The mechanical system under consideration is symmetric if its Lagrangian C{x, x) is 
invariant with respect to the action of a certain Lie group Gq on TQ, i.e., C{x,x) = 
C{gox, 9ox) for go £ Gq. In the special case when the Go-action on the tangent 
bundle TQ is the lift of a Go-action on Q (see [20]), the momentum map J : T*Q — )■ Qq, 
where T*Q is the cotangent bundle of Q and Qq is the Lie algebra of Gq, is given by 

{{Jix,V),0) = {'P,M). (3.3) 

for V := ^ the momentum conjugate to the microvelocity x- In formula (3.3), ((•,•)) 
denotes the pairing between vectors in Qq and their duals in Qq, (-, ■) is the vector-covector 
pairing on Q, and ^q(x) denotes the action of the infinitesimal generator vector field 
C,Q G X((5) defined by ^ G 0o (corresponding to the Go-action on Q) on x G Q. The 
concept of infinitesimal generator is of paramount importance in geometric mechanics 
and it is defined as the derivative of the group action. For example, the infinitesimal 
generator of the 5'0(3)-action on the initial director field no G 5*^ is 



d 
dt 



(x(^)no) = x(0)no = Uq x Hq = 2/52 (no) 



t=o 



where Uq = x(0), x(0) = I, and we omitted spatial dependence. 

In the context of Ericksen-Leslie nematodynamics, the configuration manifold is 
J-'(V, S0{3)) and the symmetry group is the subgroup Go C J-'(V, S0{3)) that leaves 
no invariant, that is, the group of planar rotations on the plane perpendicular to no. 
Upon recalhng the spatial dependence no = no(x), we identify Go = J^{T>, 5*0(2)). This 
is a typical situation in which the (gauge) 5*0(3) symmetry is broken by the 50(2) 
(gauge) group. In more generality, systems with broken symmetries involve a Lie group 
configuration manifold Q = G and a symmetry subgroup Go C G. Then, at the Lie 
algebra level, there is a subgroup inclusion i : Qq q so that the infinitesimal generator 
of the Go-action on G is induced by the infinitesimal generator of the action of G on 
itself. Upon fixing ^ G go, we have 

and on the right hand side we consider the infinitesimal generator of the right multipli- 
cation in G. Then, (^(0)^ (^) ~ X^(0 and, by denoting V = 6C/6x, the momentum 
map formula (3.3) gives 

{{j{x,v),o) = {v,m)Gix)) = {v,x^m = {{x-'vAO)) = anx-'-p) ,0) 

where the third step is due to the invariance of the pairing and the notation ((■,■)) stands 
for the pairing on g* x or Qq x go, depending on the context. The momentum map is 
therefore 

'^i.^' ^x) ^ ^X, 

and upon using the right trivialization (x, x) ^ iXi^) (with z/ = xx~^)i the above 
formula is rewritten in terms of the reduced Lagrangian as 

,f,«UrfAd;a.rfx-fxV (3.4) 
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where Ad* stands for the group coadjoint operation on G; see e.g., [20]. The latter 
formula provides Noether's conserved quantity for an arbitrary system with broken sym- 
metry. 

When the above momentum map is specialized to the case of Ericksen-Leslie ne- 
matodynamics, it is convenient to assume no = 63 so that the Lie algebra inclusion 
i : C°°(P) --)> J^{V,M.^) is i(r(x)) = (0,0, r(x))^ and the momentum map determines 
the conservation law (3.1). On the other hand, when the wryness tensor 7 is introduced, 
the choice 7o = leads to the subgroup Gq = 5*0(2) C J^{V, S0{3)), thereby leading 
to the conserved integral in (3.2) . 

4 Euler-Poincare reduction for micropolar theory 

Nematic systems are usually described by the Ericksen-Leslie theory. However, this the- 
ory limits itself to describe systems of uniaxial molecules in the absence of orientational 
defects (disclinations). As we saw, disclinations can be incorporated into the Ericksen- 
Leslie description by assuming an initial configuration 79 7^ 0- However, the resulting 
description in terms of the director n still provides no information about other possible 
molecule conformations, which can be for example biaxial. Eringen's micropolar theory 
[7] accounts for more general molecule shapes by attaching to each particle a microin- 
ertia tensor j. The resulting description has to be consistent with the Ericksen-Leslie 
theory for uniaxial molecules and this has always been an issue in liquid crystal mod- 
eling. As we shall see, the two theories are indeed consistent and, more particularly, 
Eringen's theory reduces to Ericksen-Leslie for the case of uniaxial nematic molecules in 
the absence of disclinations. 

4.1 The Lhuillier-Rey Lagrangian 

In order to insert a microinertia tensor in the dynamics, let us rewrite the primitive 
Ericksen-Leslie Lagrangian (2.5) as 

i{u,n) = - I - {u ■nf)d?i^- I F(n,Vn)d='x, (4.1) 

2 Jt> Jv 

and notice that, if we introduce the tensor parameter j(x, t) := J(I — n(x, t)n^(x, t)), 
we obtain a new Lagrangian 

£(i/,j,n) = ^ / u-jvd'^- [ F(n,Vn)d=^x, (4.2) 
^ Jv Jv 

where the microinertia tensor j(x, t) evolves according to j(x, t) = x(x, t)jo(x)x~^(x, t), 
withjo(x) := J(I— no(x)nQ (x)). Then, upon making use of the variation 5j = [Sxx~^yj]j 
the Euler-Poincare variational principle 

5 / i{u,j,n)dt = 
Jti 
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yields the equations of motion 



d 6£ _ 6i 
dt Su Su 

dtj + [j,u] = 

dfU + n X u = 0, 



+ n X — 
On 



(4.3) 



where we have introduced the notation = eihk^hk- 
the expression of i, the first equation becomes 



After some simphfications, using 



d_ 
di 



h X n . 



Then, upon allowing for more general forms of the microinertia j (which then becomes 
an independent variable), the above equations form the Lhuiller-Rey model [19] for 
continuum media. This model evidently reduces to Ericksen-Leslie when the microinertia 
tensor is given by j = J(I — nn^), which is a relation preserved by the dynamics; see 
[9]. 



4.2 Eringen's micropolar theory 

At this point, it is important to notice that, upon following the same reasoning as in the 
previous sections, one can set an initial condition of the type Vng = Hq x 7q, so that 
Vn = n X 7 at all times, provided 7 = —(Vx)x^^ + X7oX~^- This operation would 
modify the Luhiller-Rey Lagrangian by inserting the potential 7 in the set of dynamical 
variables to produce a Lagrangian 

^(i^,J,n,7) = i / ly-jud^x- / F(n, n x 7) d^x . 
^ Jv Jv 

Now, before looking into the Euler-Poincare equations for this Lagrangian, a crucial step 
consists of the observation that for j = J(I — nii"^), the Frank free energy can be written 
as a free energy functional \l/(j, 7) = F(n, n x 7) = F(n, Vn) (see [10] for the explicit 
derivation), that is (see Appendix A), 

^(J,7) = ^Tr(j7) + ^(Tr((7^)') (Tr(j) - J) - 2 Tr(j(7^)')) 

+ Tr^(j7) - ^ iv(((7jr - J^'f) , (4.4) 

where the superscript A denotes the antisymmetric part and the wryness tensor 7? is 
considered 3 real matrix, so that 77 is obtained by standard matrix multiplication 

(i.e., (7j)/fc = 'yfjak)- Notice that many terms in this expression of the free energy 
differ from those in [(3]. A particularly easy expression is available for the one-constant 
approximation (i.e., K2 = and Ku = K22 = K^^ = K): 

1 2 1 2 K K 

F(n, Vn) = -ir ||Vn|| = -fsT ||n x 7|| = — 7. ■ ^7. = — Tr(7J7) = ^(j, 7) . 
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Consequently, we obtain the micropolar Lagrangian 



whose associated variational principle gives the Euler-Poincare equations 



d 5£ S£ 
at ou ov 



6£ 



d 51 51 

+ TT-T ^ 7j X 7 — 



dtj + [j,t>] = 
. dtli + 7i y< + diV = 0. 

A short computation, using the explicit expression of i, yields the equations 



(4.5) 



d . 



d 

dx^ d^i 



l^ X 



dtj + [j, 9] = 



(4.6) 



If \l/ verifies the axiom of objectivity (i.e., \&(j, 7) = \I^(7^ ^jTZ.lZ ^^TZ) for any orthog- 
onal matrix TZ e 0(3)), then one can show (see [9]) that the first equation simplifies 

— ~ ' I X , 

dt dx^ d^,i 97" ' 

where we note that the cross product is now taken relative to the manifold indices. 
Notice that if the initial condition jo of j is such that Vjo = [jo, 7o]' then the relation 
= [j, 7] is preserved by the dynamics; this can be shown by explicitly computing 
Vj = V(xjoX^')- 

The micropolar model derived above is an extremely powerful tool in ideal nemato- 
dynamics. Indeed, this model has a twofold advantage: first, it accounts for molecules of 
different shapes and, second, it incorporates disclination dynamics that may arise from 
an initial potential 7q 7^ 0. Due to the way it is constructed, the above micropolar 
theory naturally comprises Ericksen-Leslie dynamics for uniaxial systems, which is re- 
covered upon setting Vj = [j, 7] and j = J(I — nn^) (both relations are preserved by 
the dynamics). Therefore, the above micropolar theory is the most general gauge theory 
of nematodynamics for systems composed of molecules with arbitrary shape. Indeed, 
this property follows from the fact that, if j = J{I — nn-^), then the micropolar free 
energy coincides with the Frank energy: \l/(j, 7) = -F(n, n x 7) = -F(n, Vn), so that the 
micropolar Lagrangian reduces to the Ericksen-Leslie Lagrangian. 



4.3 Remarks on biaxial nematic molecules 

This section explains how the microinertia tensor may be used to describe the dynamics 
of biaxial liquid crystals. Biaxial molecules possess two preferred directions in space, so 
that their order parameter is identified with a pair of orthogonal director fields. The 
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first dynamical description of biaxial liquid crystals appears probably in [23], where the 
equations of motion are derived by using a Hamiltonian approach and upon neglecting 
spin inertial effects. These equations were also obtained more recently in [11] by using the 
Euler-Poincare approach and upon including spin inertial effects. The treatment follows 
exactly the same steps that were described in the preceding sections, upon replacing the 
director n by the director pair (n^^^ n*^^)), so that (n*^^\ n^^^) = (xiIq^'*, xiIq^''). Then, the 
reduction procedure outlined in Section 2.2 yields the equations of motion [11] 

— — = u X — + n^'' X —- , -^— = n^'>xv (4.7) 

dt 5u 5u 5n(*) ' dt ^ ' 

where the summation convention over repeated indices is used in the first equation. The 
question of how the Lagrangian can be expressed in terms of the two director fields (and 
their gradients) is still an open question in liquid crystal dynamics. As in the uniaxial 
case, the use of the director pair (n*^^\n*^^^) for describing the biaxial molecules can 
be replaced by the definition of an appropriate tensor order parameter. In standard 
textbooks (e.g., [Jd], equation (11.9)), this tensor order parameter is identified with a 
linear combination of the following tensors 

A = I - n«n«^, B = (n« x n(2))(n« x n^^))^ - n(2)n(2)^, 

so that the total tensor order parameter is written as j = aA + /3B, where the constants 
a and /5 are thermodynamical constants. Then, it is readily seen that, the dynamical 
relation (n'^-'^^ n*^^^) = {xn^Q \ X^c)^^) implies j = xJoX'^^ where jo = o^/Kq + /3Bo is the ini- 
tial configuration involving the reference director fields (ng^'', nQ^"*). In turn, the relation 
j = xjoX~^ is precisely the same that appears in Eringen's theory and thus, if the La- 
grangian n^^\ n*^^)) in (4.7) can be expressed in terms of j as £{u, n^^\ n*^^-*) = 
then the resulting equations of motion read 



d SI SI 

— = 1/ X h 

at Sv Su 



' . Si 



Notice that the Lagrangian / involves a free energy term containing the gradient Vj. 
On the other hand, one can repeat the same steps as before, thereby introducing the 
wryness tensor (2.6) via the gauge- invariant relation Vj = [j, 7]. This last step produces 
precisely the Euler-Poincare equations (4.5), which generalize the micropolar equations 
(4.6) to arbitrary Lagrangians. In this sense, equations (4.5) were hereby shown to apply 
to either uniaxial and biaxial molecules, depending on how the microinertia tensor is 
expressed in terms of the director field(s). 

Notice that the above arguments require the underling assumption that there exists 
a Lagrangian I such that i{iy,n^^\n^'^^) = This assumption may look pretty re- 

strictive. When this assumption fails, Eringen's micropolar theory needs to be extended 
to Lagrangians of the type i{u,n^^\n^'^^) = l{iy,A, B). This step does not present any 
difficulty and was treated in [11]. 
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Another cautionary remark is also necessary here: although Eringen's microinertia 
tensor may describe the dynamics of both uniaxial and biaxial nematics, the theory 
does not account for phase transitions, e.g., from uniaxial to biaxial. Indeed, the initial 
condition iiq^'* = gives n*-^-' = x^o^'' = and B = at all times, so that if one starts 
with a uniaxial phase it is not possible to end up with a biaxial phase. 

5 Euler-Poincare approach to flowing liquid crystals 

This section reviews briefly how all previous results can be extended to flowing liquid 
crystals. In order to deal with fluid systems, it is necessary to deal with Lagrangian 
systems on the group Diff(X') of diffeomorphisms in a region V C M^, that is the Lie 
group of smooth invertible coordinate transformations on V. The next section introduces 
the topic by presenting the Euler-Poincare reduction for compressible fluids. 

5.1 Euler-Poincare reduction for ideal compressible fluids 

Two main descriptions are used in fluid dynamics: the Lagrangian and the Eulerian 
description. While the first approach describes the fluid flow by following an individual 
fluid particle, the second describes the flow as it is measured by a fixed observer in 
physical space. In the Lagrangian representation, the Lagrangian fluid coordinate is a 
(smooth and invertible) coordinate transformation r) = r7(xo,t) that maps the initial 
Lagrangian label xq to its position r] at time t. This leads to a standard variational 
principle of the type 

6 I 'c{'n,rj) dt = 

for arbitrary variations 57) that vanish at ti and ^2- The choice of the Lagrangian is 
dictated by the form of kinetic and potential energy of the fluid. For example, for 
barotropic compressible fluids, the internal (potential) energy depends only on the mass 
density, which we shall write as po(xo) in terms of Lagrangian labels. For simplicity, we 
start with pressureless compressible fluids and we shall insert the internal energy later 
on. The explicit form of the kinetic energy part of the Lagrangian is 

C'i'n^'h) = \ j^Poi^o) |'7(xo,t)|^d^xo, 

so that the corresponding Euler-Lagrange equation reads r7(xo,t) = 0, which is a free- 
motion equation. At this point, we notice that, contrary to the case of incompressible 
fluids (which are such that det{J'r]) = 1, where jT^j is the Jacobian matrix of 77), the 
relabeling symmetry is broken in the present situation: 

£(r7,r7)^£(I,77or7-i), 

where o denotes composition of maps, so that (77 o 77^^) (x) = r/(T7"^(x)). However, one 
can transfer the reasoning of the previous sections to the present case. Indeed, if we 
define the time-dependent variable 

p(x,t)= / po(xo)5(x-77(xo,t))d^Xo =: (77^po)(x, t) , 
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one can verify that 

as the following computation shows: 



^PoiV,!) = 7; [ / d^Xopo(xo)5(x-r7(xo,t)) 1(7707/ ^)(x,t)| 
Jv Jv 

= \ ^p(x) I (r/ o ^"^) (x, t) 1^ d^x 



-1^ 



= >Cp(I,77 77 

The map sending po to ry^po is known as Lagrange-to-Euler map as it maps the La- 
grangian density po(xo) to its corresponding Eulerian quantity p(x). In more geometric 
terms, one says that p is given by the push-forward of po by rj. As it happened in the 
previous sections, the Lie group action of ry on the parameter po emerges naturally in 
the reduction process and, in this case, it is precisely given by the push forward ry^po. 

Now, if we define n(x, t) := (r; o r7~^)(x, t) and i{u, p) := Cp(I,r] o r]~^), then we 
obtain the Euler-Poincare variational principle 



t2 rt2 r l^|2 



5 I i{u, p)dt = 6 I / p ^ d^x dt = 
Jti Jti Jv ^ 

subject to certain constrained variations of u and p that we now compute. First, we use 
the formula [15] 

5u = 5iriri~') = dt {i6'n)r]-') + [i5r,)ri-\ u] 

of general validity for r] belonging to an arbitrary Lie group G with Lie algebra structure 
(0, [-,•]). In the special case G = Diff(P), one has {5r])r]~^ = {5r]) o 77"^ and the 
corresponding Lie algebra is identified with the space X{'D) of vector fields, which is 
endowed with the negative Jacobi-Lie bracket 

[u, v] = (v ■ V)w — (u ■ V)v . 

Next, we compute 6p by pairing the latter with a test function ip, as follows: 

d^x5p(x)(^(x) = - / d^x(^(x) / d^xoPo(xo)V(5(x-r7(xo)) •(5r7(xo) 
V Jv Jv 



/ d^^Vfi^) / d3xopo(xo)5(x-r7(xo)) ■ (5T/ 77 i)(x) 
Jv Jv 

d^x(y9(x) div(p(x)(5T7 o r]~^)(x.)^ , 



so that 

6p = — div(p(577 o ry^"*^)) . 
In conclusion, the Euler-Poincare equations are 

-— + {u ■ V)— + divH— + Vu^- — = pV- 
01 ou ou ou ou op 

^ + div(pn) =0 
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(5.1) 



that is, after computing 51 /5u = pu and 5l/5p = 11*1^/2 and some calculations, 

du , 

— + (it ■ V)n = 
^ + div(pn) =0 

which are the equations for a pressureless compressible fluid with Eulerian velocity u. 
If we were considering (barotropic) pressure effects, it would suffice to add the internal 
energy term in the new Lagrangian 

|2 



i{u,p)= [ p\:!^d'^- [ pu{p)d?^. 

Jv ^ Jv 



Then, the Euler-Poincare equations for the velocity would be modified by the inclusion 
of a pressure term as follows: 

9u , 1 
— + {u-V)u = — Vp 

at p 

dp 

— + div(pn) = 0, 

where the pressure is given by p = p^lA'{p). Notice that the unreduced Lagrangian 
^poiVi Vl) for barotropic compressible fiuids can be reconstructed from l{u^ p) by simply 
replacing u = r} o r]~^ and p = rj^po. 

At this point, one may wonder what the correspondent of equation (3.4) is in this 
case. In the context of fiuids, the operation Ad* 5£/5z/ corresponds to the pull back 
ri*{6£/6u) and the main ingredient is the dual of the inclusion i : Xp^iV) ^ X{V). 
Here, Xp^^l^) denotes the Lie algebra of vector fields preserving po(xo)d^xo, so that 
^ G Xpg{V) <^ div(po ^) = and z(^) = Consequently, ^ G Xp^{V) if and only if there 
exists a vector potential V such that ^ = pg^curli/?, so that the Lie algebra Xp^iV) 
can be identified with the Lie algebra of vector potentials. Therefore, upon taking the 
L^-pairing of the latter relation with 'q*{6i/6u), one computes 

Then, upon considering a closed surface in V with boundary 70, one can easily apply 
Stokes' theorem to write 

0= — (h —77— ■dxo = — 0) - — ■dx = — d)w-dx, 

dt Po du dt /^(^^) p du dt 

where 7 = t7(7o) and the second step follows by changing coordinates (using the rela- 
tion p = r/^po). The above conservation relation is of paramount importance in fiuid 
mechanics and it is known under the name of Kelvin-Noether theorem. 

Remark 5.1 Notice that one can also use the L^-paring to identify the dual spaces as 
X{Vy = X{V) and Xp,{V)* = X^{V). Then we have i*{u) = poP(p^^u), where P 
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denotes the Hodge projection onto divergence free vector fields. A direct computation 
using (5.1) then shows that 

d I ..I ht 



This corresponds to Noether's theorem associated to the isotropy subgroup of po, as 
discussed in §3.3. 

5.2 Euler-Poincare reduction by stages for liquid crystals 

This section extends the previous discussion on compressible fluids to consider nematic 
liquid crystals. As we shall see, this reduction proceeds by two stages. 
For flowing liquid crystals, the simplest form of Hamilton's principle is 

5 / L(no,po)(^,'^,X,x)dt = 0, 
J to 

where r] determines a Lagrangian trajectory on the Lie group of diffeomorphisms while 
x(xo) G 5*0(3) identifies the rotational state of the point xq. The initial director field 
no(xo) and mass density po(xo) are just time-independent parameters appearing in the 
expression of the Lagrangian. By extension of the Ericksen-Leslie Lagrangian, we write 

^(no,po)(^5'7,X,x) = ^ y Po(xo) |r;(xo)|^d^xo + ^ ^ po(xo) |x(xo)no(xo) |^ d^xo 

- / po(xo) F((t7,po) o 77, (xno) o T7-\V((xno) o rf^) d=^xo . (5.2) 
Jv 

In the first stage reduction, we introduce the variables a)(xo) = x(xo)x~^(xo) and 
n(xo) = x(xo)no(xo). Then, the Lagrangian possesses the following symmetry 

However, this second Lagrangian also possesses another type of symmetry that is related 
to the relabeling properties. Upon defining the variables 



u := uj o T] ^ , n:=nor7^ 



this symmetry is 



where the action of the diffeomorphism group again emerges in the definition of the 
variables 1/ and n and 

i{u, p,u,n) = - f p|u|^d^x+— f p|i/xn|^d^x— f p F{p~^,n,Vn) d^x. 

At this point, upon denoting w = 6r] o rj^^ and ^ = {Sx)x~^ ° the variational 
principle 

/•t2 

6 / i{u, p, i^, n) dt = 

Jti 
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for the constrained variations (see [9] for the derivation of these formulas) 



5u = dfW + [w, u] 
6p = — div(piu) 

6iy = dt^ — w ■ Vv + u ■ V$, — v x ^ 
6n = —w ■ Vn + ^ x n 



yields the Euler-Poincare equations 



6u 



Vw'^ ■ = -Vu ■ - Vn^ ■ ^ 



61 
6^ 



5l_ 



pV- 
6p 



d_5i_ 

dp 
{ dt 



diiu 



51 



V X 



6I_ 

dn 

'dt 



n X 



6I_ 
6^ 



+ [u ■ V)n = 1/ X n . 



(5.3) 

Taking the variational derivatives of the Lagrangian I given above yields the Ericksen- 
Leslie equations of nematic liquid crystals 



P 



J 



du 
'dt 



phi 

Dt^ 



n ■ h + Jn 



-Vp — di p Vn 



dF 
dn,i 



phi 

Dt^ 



n 



0, 



where we have defined the pressure p := —dF/dp~^ and used the material derivative 
notation D/ Dt := dt + u - V . 

At this point, one can repeat all the steps that were followed in the fiowless case 
and insert an initial condition Vng = ng x 7. Then, the relation Vn = n x 7 is again 
invariant, provided the diffeomorphism ry acts on it by push- forward. An analogous 
result also holds for the introduction of the microinertia tensor j. At the end of the 
various steps of the reduction, one obtains a Lagrangian £(■«, i/, p, j, 7) that reads 

^{u,u,p,j,-f) = 1- I p|ii|^d^x+^ / pu-jud^^- I p^(p"\j,7)d^x, 

where the free energy is such that \E'(p~^,j, 7) = F(p~^,n, Vn) when j = J(I — nn^). 
In turn, upon making use of the index summation convention, this produces the system 



d 6i 
dt 5u 



d 5i 
di'6u 
dp 



[u ■ V — + div u — + S/v/ ■ — 
ou du du 



-S7W- — + pV— 



op ojhk hi dx^ 



V X 



dt 



+ di{u'p) = 0, 



9tJ + [j, i>] = 0. 



d 51 5£ 
dx' 57, 57, 



(5.4) 



+ 7i X + diU = . 
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whose first two equations become 



when \l/ verifies the axiom of objectivity, see [9] for details. 

Notice that, although here we presented the equations for compressible fiows, their 
incompressible version is easily derived by setting p = 1 and by adding the pressure 
force — Vp on the right hand side of the velocity equation, which then becomes 

The last step, i.e., adding the pressure force term, is quite common in fiuid mechanics and 
it amounts to projecting the whole equation to its divergence-less part. Alternatively, the 
pressure p can be inserted as a Lagrange multiplier in Hamilton's principle to constrain 
the relation p = 1; see [ ] for explicit calculations. 



5.3 Conserved quantities and helicity of uniaxial nematics 

Conserved quantities also appear for nematic liquid crystals, at least in the case of 
uniaxial molecules. For example, it is a simple verification to see that equations (5.3) 
produce the following dynamics of the momentum map n ■ 6i/6u: 

d f6e \ ^. f f6£ , , ^ 
n + div Ut — ■ n = 



dt \Su J \ ysu 

so that the corresponding integral is constant. 

On the other hand, one may ask whether more general dynamical invariants may 
appear in liquid crystal dynamics. For example, it would be natural to wonder if the 
celebrated helicity invariant of Euler's fiuid equation may have a correspondent in liquid 
crystals. A positive answer was given in [11, 12], where it was shown that 



d 
dt 



(p - \ — ■ (Vn X n) ■ dx = — (tiu-Ju- (Vn x n) • dx = 



for any loop 7 moving with the fiuid velocity u. Such a Kelvin-Noether circulation 
theorem is produced by the fiuid relabeling symmetry underlying liquid crystal dynamics 
[15]. Fluid models do not always possess such conservation laws, so that the right hand 
side is often expressed as a non-zero circulation integral (this is the case, for example, of 
magnetohydrodynamics in plasma physics). However, the case of liquid crystals benefits 
from a general theorem in geometric mechanics that is due to [17] (see Proposition 
2.2 therein). This result is based on the momentum map — (Vn x n) ■ Si/6u that is 
produced by the action of the diffeomorphisms on the Poisson manifold so (3)* x S"^, with 
coordinates (cr, n). Then, the total circulation quantity 6i/6u— (Vnxn) -Si/ 6u identifies 
the momentum map emerging from the diffeomorphism action on x so(3)* x S^, 
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which in turn generates the Kelvin-Noether theorem. See [11, 12] for further details. 
At this point, taking the curl of the total momentum map gives the expression of the 
hydro dynamic vorticity 

fl = curl(w — Jv ■ (Vn x n)) 
and one verifies that the hydrodynamic helicity is constant 

71= / [u — Jv ■ (Vn X n)) ■ curl(n — Jv ■ (Vn x n)) d^x . 

Jv 

The latter relation does not come as a surprise in geometric mechanics. Indeed, it is 
simply given by the usual expression of fluid helicity upon shifting the fluid momentum 
by the quantity —Jpv ■ (Vn x n). In turn, as noted earlier, the latter emerges from a 
momentum map; this type of shift by a momentum map agrees with the general theory 
developed in [17] (see Corollary 2.3 therein). It is not the purpose of this paper to 
present how this more involved theory works; however, it is still worth emphasizing that 
the geometric setting underlying nematodynamics allows to find conserved quantities and 
circulation theorems by simple application of general principles in geometric mechanics. 



6 Remarks on dissipative dynamics 

The dissipative case can be easily handled by using a Rayleigh dissipation function 
i?no(X)X)) ^ (see, e.g., [20, §7.8]), having the same invariance prop- 

erties as the Lagrangian >Cno(X;X) and thus producing a reduced function r{v,n) = 
r(xx~^,Xno) = -Rno(X;X)- The associated variational principle 



yields, by reduction, the variational principle 

rt2 



rt2 rt2 

6 £{v,n)dt= (jy, n) ■ dt, 

Jh Jh ^v 



(6.1) 



for variations 6n = (5x)x^^n and 6v = dt {{Sx)x^^) + [{^x)x^^:'^]- We thus obtain 
the reduced Euler-Poincare equations with dissipation (these are the reduced Lagrange- 
d'Alembert equations when the external force is dissipative and the configuration space 
is a Lie group) 

d 6£ 6£ 6i 6r , 

— 7 — h — xi>' = nx — (i/,n) 

Ot OV OV on QV (r. o\ 

^ (6-2) 
an 

— + n X 1/ = 0. 

at 

Dissipation in Eringen micropolar theory can be handled similarly, by considering 
a Rayleigh dissipation function i^^g ,y|j)(x, x), invariant under the isotropy subgroup of 
(jo,7o)) thus yielding a function 

r{y,hl) = r (xX"\xJoX"\-(Vx)x"^ + X7oX~^) = ^(jo,7o)(^' 
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Applying the same variational principle as above, we get the following modification of 
equation (4.5): 



5^ 

X — + 



d 5i 
dt3 + [3, P] = 

dtii + ji y< + 



J, 



6i 



0. 



+ 



d 5i 



5t 5r , 



The case of flowing liquid crystals can be handled similarly. The balance of momentum 
and the balance of moment of momentum equations in (5.4) are modified by the addition 
of the terms 

-^(u, i^,p, j,7) and - — (u, i/, p, j, 7), 



6u 



6iy 



respectively. 

Remark 6.1 (Dissipation function) Notice that the above equations involve an ar- 
bitrary dissipation function r. For example, this is the approach followed in fundamental 
treatments of liquid crystal dynamics, e.g., [1]. However, it would be interesting to spe- 
cialize the above relations to a particular expression of r and study how one can pass from 
one representation to the other, as it was done for F(n, Vn) and \l/(j, 7) in Appendix 
A. We leave this interesting question as an open subject for future work. 
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A The free energy 

We reproduce here the computation in [lit] that shows how all terms in the Frank energy 
F(n, Vn) := K2 (n ■ curln) +]-Ku (div n)^ +^ir22 (n • curln)2+^K33 ||n x curln|p, 

chirality splay twist bend 

can be rewritten in terms of the variables j = J(I — n n) and 7, the latter being 
introduced through the invariant relation Vn = n x 7. Thus, the explicit expression for 
the micropolar free energy \l/(j, 7) of nematic media given in (4.4) will be written after 
computing the micropolar expression for each term in the Frank energy. Further details 
are given in 

Twist. Using n (g) n = I — j/ J, we have 

n ■ V X n = — n ■ 7(n) + ||n|p Tr(7) = — Tr (77) = — Tr {jj^) ■ 

where 7(n) = -iiaUa, with a being the so (3) ~ M'^-index, and 7'^ denotes the skew part 
of 7, i.e., 7'^ = (7 — 7"^) /2, where we see 7 as a 3 x 3 matrix with components 7^^. 
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Splay. We introduce the vector % = eabdaci defined by the condition 7-u = Tr(u x 7), 



for all u G 
compute 



(u X -iX- We 



where u x 7 is the matrix with components (u x 7)^^ 
(div n)^ = (7 ■ n) (7 ■ n) = 7 ■ (n (g) 11)7 

= 2 (Tr(j)/J - 1) Tr ((7^)^) - ^ Tr {j{^^f) , 

where 7"^ denotes the skew part of 7, i.e., 7^ = (7 — 7-^)/2 and where we used the 
equality 7 = —27"^. The latter can be shown by noting that we have the equalities 
Tr(^u) = -27 ■ u = -2 Tr(7u) for all u G Ml 

Bend. For all u G M^, we have 

(n X (V X n)) ■ u = — Vnn ■ u = — (n x 7(11)) ■ u = — (u x n) ■ 7(11) 
= — un • 7(n) = — Tr((un)"^7n) 
= Tr ((n (g) n){i7) = Tr ((I - j / J)u-i) 

= Tr(uX7)-^Tr(ux (7j)) 



so we get 



and therefore 



7-u - -13 ■ u. 



n X (V X n) = 7 — -7j 
u 



|n X (V X n) 



-jl3 -7 



-2Tr( (-(7jy-7^ 



Summing all the terms, we obtain that the Frank free energy is indeed given by the 
expression for \E' provided in (4.4). 
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